BIRTH AND DEATH IN DISCRETE MORSE THEORY 
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Abstract. Suppose M is a finite simplicial complex and that for = to, t\, . . . , t r = 
1 we have a discrete Morse function F ti : M — > R. In this paper, we study the births 
and deaths of critical cells for the functions F ti and present an algorithm for pairing 
the cells that occur in adjacent slices. We first study the case where the trian- 
gulation of M is the same for each f,-, and then generalize to the case where the 
triangulations may differ. This has potential applications in data imaging, where 
one has function values at a sample of points in some region in space at several 
different times or at different levels in an object. 



I. Introduction 

The purpose of this paper is to study the discrete analogue of the following phe- 
nomenon in classical smooth Morse theory. Suppose that N is a smooth manifold and 
that we have a family of functions / : N x / — » R such that the various f t : N — > R 
are generically Morse; that is, for almost all t, the function f t has only nondegener- 
ate critical points. Then as t varies, the critical points of the f t move around in N. 
Sometimes a critical point is "born"; that is, a new critical point appears at some 
time to- At other times, critical points "die". Generically, critical points are born 
and die in pairs. Such events are isolated since the critical points of a Morse function 
are separated; we call such points in N x / birth-death points. 

Consider the following simple (but non- generic) example. Let N = R and consider 
the family f : N x [-1, 1] -> R defined by 

4 3 

f(x,t) = — + j +tx 2 . 

Each ft is a Morse function except for t = where x = is a degenerate critical 
point, and t = 1/8 where a degenerate critical point occurs at x = —1/2. Figure 
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Figure 1.1. The curves x 4 /4 + x 3 /3 + tx 2 for t = -.9, 0, 3/32, 1/8, .9. 




Figure 1.2. The bifurcation diagram associated to x A /4 + x 3 /3 + tx 2 . 

1.1 shows the graphs of f t (x) for a few values of t. Note that for t < 0, there are 
three critical points of f t , two of index and one of index 1. At t — 0, the critical 
point x = (—1 + y/l — 8t)/2 that existed for t < dies, and we see the emergence 
of the degenerate critical point x = 0. Then for < t < 1/8, we again see three 
critical points, but now two of them occur at negative values of x. At t — 1/8, 
the two negative critical points come together to form a degenerate critical point at 
x = —1/2. Finally, for t > 1/8, we have only the critical point x = 0, which is critical 
for all t E [-1,1]. 

We can trace the critical points as t varies via a bifurcation diagram as in Figure 
1.2. Note that in this simple example, the diagram is merely the graph of f x (x, t) = 
with x on the vertical axis and t on the horizontal. The two components of the curve 
intersect at t — 0. 

We are concerned here with the following question. Suppose we have a simplicial 
complex M and for each = t ,ti, . . . ,t r = 1 we have a discrete Morse function 
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F ti : M — > R (the basics of discrete Morse theory are given in Section [2]). We allow 
the possibility that we have chosen a different triangulation of M for each i — 0, . . . , r. 
We wish to study how the critical cells move in M as % varies and construct bifurcation 
diagrams similar to that in Figure 1.2. To do this, we need to extend the discrete 
Morse functions F ti to such a function on the complex Mxl and develop an algorithm 
to pair critical cells in each slice with those in the previous one. We carry this out 
first in the case where the triangulation of M is the same for each F tl (Section [3]) and 
then extend this to the general case (Section E]). 

In Section [S], we show the results of running our algorithms on CT scans of the 
head, where the function in question is density. Minima correspond to cavities in the 
interior of the head. Also, in a series of renal scintigrams we trace the maxima of the 
intensity of radiation which correspond to the locations of the kidneys and bladder 
within the image. 

1.1. Acknowledgements. The results in this paper were obtained at the Workshop 
in Computational Topology held in conjunction with the 20th Summer Conference 
on Topology and its Applications at Denison University, July 10-14, 2005. We thank 
Tom Peters for organizing the workshop, thereby providing an excuse for us to meet 
again. We also thank Denison for its hospitality. 

2. Discrete Morse Theory 

In this section we review the basics of Forman's discrete Morse theory |F2J. Let 
M be a simplicial complex. A typical p-simplex will be denoted by and if a is 
a codimension-one face of r, we write a < r. In the following definition, we shall 
abuse notation slightly by speaking of functions from M to R when we really mean 
functions from the set of simplices of M to R. No confusion will result. 

Definition 2.1. A function / : M — > R is a discrete Morse function if for every 
o<vp) (= M, the following two conditions hold: 

(1) #{^) > a\f{(3) < /(a)} < 1, 

(2) #{ 7 (^ 1 ) < a\f(j) > /(a)} < 1. 

Essentially, then, discrete Morse functions are functions on M that increase with 
the dimension of the simplices; that is, the values on all but at most one face of a 
must be smaller than the value on a itself. Discrete Morse functions exist. Indeed, 
the simplest (and most trivial) example is the following. If a is a simplex of M, we 
define / : M — > R by 

f(a) = dim a. 

Note that the conditions (1) and (2) in the definition of a discrete Morse function 
are exclusive. That is, if one of the sets has cardinality one, the other is empty. To 
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see this, assume that E M has a coface t^ p+1 ^ such that /(r) < f(o~) and a face 
v^' 1 ' such that f(v) > /(cr). Let tr'W be a different face of r which also contains 
t>. Then the value /(cr') must be bigger than /(f), since v already has a coface with 
smaller value, and similarly the value of r must be bigger than /(cr'), so 

/(r) < /(a) < /(f) < /(a') < /(r) 

which is not possible. 

Smooth Morse theory depends heavily on the notion of critical point. The discrete 
theory does as well. 

Definition 2.2. Let / : M — > R be a discrete Morse function. A simplex is 
critical if the following two conditions hold: 

(1) > a\f{(3) < /(a)} = 0, 

(2) #{7^ < a\f(j) > /(a)} = 0. 

A simplex that is not critical is called regular. 

For example, if / : M — > R is given by /(cr) = dim cr, then every cell is critical. 

Theorem 2.3 QF2J, Theorem 2.5). Suppose M is a simplicial complex with a discrete 
Morse function. Then M is homotopy equivalent to a CW-complex with exactly one 
cell of dimension p for each critical simplex of dimension p. 

Constructing discrete Morse functions is rather involved, especially if the complex 
has many simplices. Recall that in smooth Morse theory, one often works with the 
gradient vector field of a Morse function. This encapsulates the qualitative behavior 
of the function so that one need not even know the values of the function itself. In 
this paper, we shall in fact construct the gradient of some discrete Morse function. 

The fundamental observation is this: regular simplices occur in pairs. For an 
arbitrary simplicial complex with a discrete Morse function /, we draw arrows to 
represent a vector field as follows. If aS^ is a regular simplex with p( p+1 ^ > a satisfying 
f{/3) < f{ct)i then we draw an arrow from a to (3. Figure 2.1 shows an example of 
a vector field on the torus (identify the top edge with the bottom, and the left edge 
with the right). 

It is easy to see that every simplex a satisfies exactly one of the following: 

(1) a is the tail of exactly one arrow; 

(2) a is the head of exactly one arrow; 

(3) a is neither the head nor the tail of an arrow. 

A simplex is critical if and only if it satisfies condition (3) above. These arrows can 
be thought of as the gradient vector field of the Morse function. A better point of 
view is the following. A discrete vector field V can be thought of as a collection 
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Figure 2.1. A discrete vector field on the torus 

pairs {aS p \ p( p+1 ^} of simplices, where a pair a, [3 is in V if and only if a < /3 and 
f{0) </(«)• 

Definition 2.4. A discrete vector field V on M is a collection of pairs {a^,/3^ +1 ^} 
of simplices of M with a < /3 such that each simplex is in at most one pair of V. We 
often use the notation V(a) = f3 to indicate that the pair {a,/3} G V. 

Of course, one would like to know if a discrete vector field V is the gradient field 
of some discrete Morse function on M. The following definition is crucial. 

Definition 2.5. Let V be a discrete vector field on M. A V-path is a sequence of 
simplices 

such that for each i = 0, . . . ,r, {a r ,/3 r } G V and /3j > a^+i ^ «j. Such a path is a 
non-trivial closed path if r > and «o = «r+i- 

Forman proved the following two theorems. 

Theorem 2.6 (|F2J, Theorem 3.4). Suppose V is the gradient vector field of a discrete 
Morse function f . Then a sequence of simplices is a V-path if and only if di < Pi > 
aj+i for i — 0, . . . , r, and 

/(« ) > /(A)) > fM > f(Pl) > ■ > /(ft) > /(Or+i). 

In particular, if is a gradient vector field, then there are no nontrivial closed 
V^-paths. The converse is true as well. 
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Theorem 2.7 QF2J, Theorem 3.5). A discrete vector field V is the gradient vector field 
of a discrete Morse function if and only if there are no nontrivial closed V -paths. 

The reader is invited to check that the discrete vector field shown in Figure 2.1 is 
the gradient of a discrete Morse function on the torus. 

A convenient combinatorial description of vector fields may be given in terms of 
the Hasse diagram of M. This is the directed graph whose vertices are the simplices 
of M, and whose edges are given by the face relations in M (i.e., there is an edge 
from (3 to a if and only if a is a codimension-one face of (3). Given a vector field V, 
modify the Hasse diagram in the following manner. If {a < /3} G V, then reverse the 
orientation of the edge between a and (3. A V^-path is then a directed path in the 
modified graph. 

Theorem 2.8 (|F2J, Theorem 6.2). There are no nontrivial closed V -paths if and only 
if there are no nontrivial closed directed paths in the modified Hasse diagram. 

3. The Algorithm, Part I 

As a special case, in this section we assume the following. Let M be a finite 
simplicial complex and suppose that for each = to, t±, . . . , t r = 1, we have a discrete 
Morse function F t . : M — > R defined using the same triangulation of M for each %. 
In applications, we would probably have the values of each F ti only on a collection of 
points in M, but we can always extend this to a discrete Morse function on all of M 
via the algorithm we presented in [KKMj . We therefore assume that this has already 
been accomplished. 

We shall call the complex M x {tj} c M x I a slice. Our algorithm constructs 
a discrete vector field on M x / by successively adding each slice and extending the 
existing field on M x [0, tj_i] U M x {ti}. In this case, where the triangulation of each 
slice is the same, we shall not end up with a triangulation of M x /, but rather a 
regular cell decomposition. The discrete Morse theory is the same in this case; see 
[FT]. 

The algorithm is quite straightforward. In the space M x [0, 1], we have a trian- 
gulation of each slice M x {ti}. Create a regular cell decomposition of M x [0, 1] by 
adding the cells a x [U-\, U] to the slices. Suppose we have defined the vector field V 
on M x {ti-i}, and we know the vector field associated to F ti . We denote the vector 
field on the slice by Vi. We therefore have pairings of cells in the slice M x {U}, with 
some cells being critical for the discrete Morse function F ti . These critical cells are 
not paired with a cell of higher dimension by Vi. The procedure is then as follows. 

• For each critical cell o~GMx{tj}, define V{o~) = o~X [ti-i, ti] ; that 
is, pair a with the cell a x [ti-i,ti\. 
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Figure 3.1. The vector field on S 1 x I 
• For each regular cell r G M x {U} , define V[r x = Viij) x 

Theorem 3.1. The resulting vector field V on M x [0,t,] is the gradient field of a 
discrete Morse function. Moreover, the procedure does not create any new critical 
cells. 

Proof. This is more or less clear. Since the vector field Vi has no nontrivial closed 
paths, the extension of V to M x [0,U] cannot either since the arrows always point 
"backward" in the t direction. That is, there could be no closed path of simplices in 
M x since the field V is basically Vi x id there, and since we may assume 

inductively that there are no closed paths in M x [0, and no path entering from 
M x can return, we are finished. As for critical cells, those cells that are 

critical for F ti are now paired in the extended field V, and no new critical cells occur 
in M x [ti-i,ti\. Note that this is in keeping with Theorem 12. 31 □ 

Figure 3.1 shows an example in the case M = S 1 . There are eight slices in the 
vertical direction. In each slice the vector field Vi is represented by black arrows and 
the critical cells are shown in gray. The extension of the field to each M x ti\ is 
also shown with black arrows. 

4. Birth and Death 

The procedure described in Section [3] gives us a mechanism to follow the births 
and deaths of critical cells in the various slices. Observe that as we extend the vector 
field V from M x [0,tj_i] to M x [0,U], we create a V^-path from the cells that were 
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critical in the slice M x {tj} to a K-path in the slice M x {tj_i} which must lead to 
some critical cell of the same dimension. This allows us to track births and deaths of 
critical cells in M as t increases, similar to the example presented in the introduction. 
See Figure 4.1 for an example in the case M = S 1 (here, the top and bottom vertices 
of each slice are identified). 

A bit of explanation is in order. We have labelled the vertices 1 through 8 and the 
edges a through h. Consider first the critical vertices. Given a critical vertex v in 
slice i, the extended vector field pairs it with v x [ti_i,ti\. The other vertex of this 
edge is then a part of a path for the vector field V^_i, and we follow it to a critical 
vertex. This is how we pair v with a critical vertex in slice i — 1. Note that there is 
a birth in slice t\ and that persists until slice £3, when the cell dies. 

The critical edges are a bit more complicated. A \^-path beginning with a critical 
edge e then passes to the cell e x [ti-i,ti\. From there, we have two directions in 
which we may proceed. In the bifurcation diagram, if a path starting with e contains 
a cell which has a critical edge in slice i — 1 as a face, we join e with that edge. It is 
possible that neither direction passes by a critical edge; this corresponds to a birth. 
If both directions pass by a critical edge, we get a death. 

Note that we may also run our procedure in backwards time as well. In forward 
time, each critical vertex at time t is connected to a unique critical vertex in time 
t — 1. However, it is possible that several in time t are connected to the same one 
in time t — 1. When running the algorithm backwards, though, each critical vertex 
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Figure 4.2. The vector field on S 1 x / in backwards time 

in time t — 1 is connected to a unique critical vertex in time t. This leads to the 
following definition. 

Definition 4.1. Two critical vertices are strongly connected if they are connected 
by a path in both directions. For critical edges and higher dimensional critical cells, 
which might be connected to several other such cells in the previous slice, and in fact 
might have several paths in space-time between critical cells, we require a single path 
in space-time in each direction. 

Figure 4.2 shows the result of running our algorithm backwards on the data shown 
in Figure 3.1. The resulting bifurcation diagrams are presented in Figure 4.3. 

A quick study of the diagrams in Figure 4.3 reveals the strong connectivity diagrams 
shown in Figure 4.4. 

5. The Algorithm, Part II 

We now consider the general case where we allow different triangulations of M in 
each slice. Indeed, this case is important for applications where the sample of points 
taken from an object may change over time, thereby yielding different triangulations. 
Since these samples form the basis of our algorithm for generating discrete Morse 
functions |KKMj . we must consider this possibility. 

For each i = 0,1, ... ,r, denote by Mj the complex M with its associated trian- 
gulation, and suppose for = t ,ti, . . . ,t r = 1 we have a discrete Morse function 
F ti : Mi — > R. Let Vi be the vector field on the slice = M x {ti} associated to F ti . 
We now proceed to define a vector field on M x I extending the various VJ. 
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ti t-2 t3 t4 ts tg t7 ti t2 t3 t4 ts t6 t-j 

Figure 4.3. The births and deaths of minima and maxima in back- 
wards time 



strongly connected minima strongly connected maxima 




tl t2 t3 t4 ts tg ti ti t-2 t3 t4 ts t-c u tj 



Figure 4.4. The strong connectivity diagrams 

Create a triangulation of Mx[0,ij by extending that of Mx[0,tj_i] 
to M x [0, ti-i] U Mi . Note that this is not unique. 
For each critical p-cell a G M i} consider the collection of (p + 
l)-cells Vi * a, . . . ,vg * a , V{ G Mj_ ls having a as a face. Find the 
cell r = f * a, where v is closest to a critical p-cell in Mj_ ls 
and set V^c) = r . 




• Extend the vector field to the remaining cells so that no closed 
V-paths result; this is possible by Lemma 4.2 of IIF1II . 

Some remarks are in order. Note that this procedure is rather inefficient. First, the 
triangulation we choose is not unique. This is not a serious problem, but it does leave 
some ambiguity. Our choice for V(a) for critical cells a G Mj is based on the idea 
that minimizing the length of the path out of the cell should lead to the critical cells 
in Mi_i "closest" to a in M. Finally, observe that the extension to all of M x [0, tj\ 
may have critical cells in M x This is not serious, and we can cancel as 

many of them as possible via the cancellation algorithm presented in |KKM] . Note 
that, presumably, we should be able to cancel all of the new critical cells that occur 
since M x [0,ti] collapses onto M x [0, ij-i]. However, we do not claim to have an 
algorithm which guarantees this. 

Figure 5.1 shows an example with M — S 1 . The critical cells in each slice are 
shown in gray. Note that our extended field has no additional critical cells. 

6. Birth and Death, Part II 

The procedure described in Section H] works equally well in the general case of 
different triangulations of M in each slice. The only modification is that there is no 
consistent labelling of the simplices that we may use on the axes of the bifurcation 
diagrams. Still, in each slice, we may enumerate the critical cells and keep track of 
births and deaths as before. Figure 6.1 shows the diagrams associated to the vector 
field in Figure 5.1. 
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minima maxima 



ti tz t3 ti ts tg ti t 2 t 3 t 4 t 5 t 6 

Figure 6.1. Minima and maxima for different triangulations of S l 

It is possible to speak of strong connectivity for critical cells in this context as well. 
The interested reader is invited to construct the associated diagrams for the example 
above. 

7. Relation with smooth Morse theory 

Suppose /: M — > R and g: M — ■> R are Morse functions and g is close to /. If 
p G M is a critical point of / with index i there will be an index i critical point 
q of g close to p. Let D?(p), D 9 (q), A^(p), A 9 (q) be the descending and ascending 
manifolds of p and q, i.e., if (f>t(x) is the flow of a gradient-like vector field for / and 
ipt{x) is the flow of a gradient-like vector field for g then 

D^(p) = {x E M | lim (j>t(x) = p] 
D 9 (q) = {x G M | lim ip t (x) = q] 

t^oo 

A (p) = {x G M | lim =p} 

oo 

A 9 (q) = {x G M | lim ^(x) = g} 

t— oo 

We know that D^{p) and ^4^(p) intersect transversely in the single point p. Likewise, 
D 9 (q) and A 9 (q) intersect transversely in the single point q. Since / and g are close we 
know that D*{p) and D 9 (q) are closej. Hence D^(p) and A 9 (q) intersect transversely 



1 or to be more precise, they are close in some neighborhood of q. If x € (p) then 4>t{x) € (p) 
also, but you might not expect (j>t{%) to be near to D 9 (q) for large negative t. 
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in a single point r and D 9 (q) and A* (p) intersect transversely in a single point s. 
This means that there is a unique path from p to q which first travels along the flow 
and then along the flow ip. (The first part of the path is 0t(r) for alH > and the 
second part is ^(r) for all t < 0.) There is likewise a unique path from q to p which 
first travels along the flow if) and then along the flow 0. Our discrete scheme mimics 
this, using V-paths instead of the flows. 

It is natural to ask to what extent the birth-death diagrams our algorithm produces 
agree with those in the smooth case. That is, if we have a family of smooth maps 
G : iV x / — > R which are generically Morse, and after choosing a triangulation of N 
and restricting the maps G(—,t) to a finite collection of slices we run our algorithm, 
how is the resulting birth-death diagram related to that for Gl To begin the study 
of this we first ask the following question. Given a Morse function F on a manifold 
M, does there exist a triangulation of M and a discrete Morse function on that 
triangulation whose critical cells are in one-to-one correspondence with the critical 
points of Fl This is answered by the following theorem of Gallais [G] . 

Theorem 7.1. ([G] Theorem 3.1]) Let M be a smooth Riemannian manifold and 
let F : M —>■ R be a Morse function. Then there is a C 1 -triangulation of M and a 
discrete vector field V such that the critical cells ofV are in one-to-one correspondence 
with the critical points of F. Moreover, the critical points lie in the interiors of the 
corresponding critical cells. 

Given a compact smooth manifold M and a Morse function / : M — > M one might 
expect that for any triangulation K of M with sufficiently small mesh there is a one- 
to-one index preserving correspondence between the critical simplices and the critical 
points of /. The following example shows that this is far from true. 

Consider a compact surface M and a smooth map / on M which, in local co- 
ordinates around some point p G M, is given by / : M. 2 —>■ K given by f(x,y) = 
(y — 10x)(y — llx). Then / has a single critical point, the origin, and it has index 
1. Now for any < e < 1 take a triangulation of IR 2 with vertices on the lattice 
{(em, en) | m, n G Z}. Then all vertices in the link of the origin have values greater 
that so the origin is a critical vertex, although / does not have a critical point of 
index near the origin. In fact, the closest point on the lattice with negative value 
of / is (2e,21e). 

This example shows that for a fixed triangulation additional critical simplices not 
corresponding to critical points of / may be unavoidable, no matter how small the 
mesh. Moreover, the measure of the set of small mesh triangulations with extra 
critical simplices is large. 

Theorem 17.11 allows us to make some remarks about the birth-death diagrams 
generated by the algorithms above. 
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Theorem 7.2. Suppose M is a smooth Riemannian manifold, and let G : M x I — > R 

6e a family of smooth functions such that the G(—,t) are generically Morse. Choose a 
finite collection of values = to < t x < ■ • • < t r = 1 with each G ti '■— G(—, tj) Morse. 
Then if the ti are sufficiently close, the birth-death diagram produced by algorithm #2 
is a piecewise linear approximation of the birth-death diagram of G. 

Proof. By Theorem 17.11 for each i = 1, . . . , r there is a triangulation of M and a 
discrete Morse function on Ki with critical cells in one-to-one correspondence with 
the critical points of G^. Now consider the result of applying algorithm #2 to this 
collection of discrete Morse functions. If the points ti are sufficiently close together, 
then the critical points of Gt i+1 will be close to those of G^', however, there may be 
births or deaths in the interval (tj, tj+i), and so there is not necessarily a one-to-one 
correspondence between the two sets of critical points. Moreover, the triangulations 
Ki and K i+ i will be similar in the sense that there is a flow in the space of C 1 - 
triangulations on M taking Ki to Ki + i, and since the ti are close, the triangulations 
are close in this space. Thus, the resulting discrete Morse functions fa and f\ + \ 
will have their critical cells close together. Algorithm j^l will then pair the correct 
simplices, and the resulting diagram is therefore a piecewise linear approximation of 
the birth-death diagram for G. □ 

Note, however, that if the points ti are too far apart, the birth-death diagram may 
pair the "wrong" critical cells. For example, consider the torus T, which is the surface 
of revolution obtained by revolving the circle x 2 + (y — 3/4) 2 = 1/16 about the x- 
axis, with the height function given by projection to the z-axis so that the absolute 
minimum is z = — 1 and the maximum is z = 1. Now consider the family of functions 
G : T x [0, 27r/3] — > R obtained by taking the height function given by projection to a 
line in the yz-pl&ne obtained by rotating the z-axis counterclockwise about the origin 
through an angle 9 for < 9 < 2tt/3. The index 1 critical point p = (0, 0, 1/2) at the 
top of the hole will then move along the curve y 2 + z 2 = 1/4 tracing the interior of 
the hole, ending at the point q = (0, — v3/4, —1/4) which is the critical point at the 
"top" of the hole for G 27T /3. Let r = (0,0, —1/2) be the critical point at the bottom 
of the hole for Go; the point — q will be at the bottom of the hole for 6^/3. In the 
birth-death diagram for G, there will be a smooth curve joining p to q, and a smooth 
curve joining r to —q. But, if we consider the discrete version of this problem with 
to = and ti = 27r/3, then we will find that the critical edge corresponding to p will 
be joined with the edge for —q, and the edge for r will be joined with that for q. 
With care, we can even use the same triangulation for t and t\ and use algorithm 
#1, obtaining the same result. 

In practice, however, this is unlikely. After all, discrete approximations are usually 
chosen to be sufficiently dense, and so with care the above situation can be avoided. 
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8. Applications 

Birth-death diagrams produced by this algorithm can be used for tracing features 
through a sequence of images. In this section we describe two such applications. Our 
purpose is not 3D reconstruction, though, since a number of reconstruction algorithms 
(for example |NCMO] ) and even patents, e.g. 

http:/ /www.freepatentsonline. com/7079132, html 

and 

http:/ /www. freepatentsonline.com/6606091. html 

exist. 

The algorithm was tested on two sequences of images. The first sequence is obtained 
from a CT scan of the head where the algorithm is used to trace minima of the density 
function which correspond to cavities and tunnels. 

The number of local minima of the density function in one image ranges from 200 
to 300, and most of these correspond to small fluctuations in the density and to noise. 
In order to eliminate these, as well as local minima with higher function values which 
correspond to softer tissue, we use persistence. Persistence is incorporated in the 
algorithm for generating a discrete vector field and enables cancelling pairs of critical 
simplices that are connected by just one gradient path and have values that differ by 
less than the chosen persistence level p. By setting the persistence level to a suitable 
value, pairs of critical points which correspond to a small fluctuation in the density 
function, as well as local minima which represent softer tissue are cancelled, while 
local minima with low values of the density function which correspond to the cavities 
are preserved. 

Figure I8TT1 shows the birth and death diagram of the minimal vertices in a sequence 
of 12 images obtained by running the algorithm with persistence set to p = 30 (where 
the values range from to 255). The forward connections correspond to thin lines, 
while the strong connections correspond to the thicker grey lines. 

Figure 18.21 shows the sequence of images in which the cavities connected by strong 
connections are traced. Points which are strongly connected are labeled by the same 
letter. Altogether 12 strong connections have been found. As can be seen, each strong 
connection represents a cavity, and most cavities are detected. Figure 18.31 shows the 
connections among the local minima, simulating the cavities themselves. 

We have also tested the algorithm on a series of scintigrams obtained in renal 
scintigraphy. Images of the abdomen are taken with a gamma camera at fixed time 
intervals after a radioactive agent has been injected. The purpose is to follow the flow 
of the agent through the kidneys in order to evaluate their function. The algorithm 
was used to trace points of maximal concentration of the agent which correspond to 
locations within the kidneys and the bladder. 
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Figure 8.1. Births and deaths of minima in a sequence of head slices 

In contrast to CT scans, where noise is present in the form of a relatively small 
error in the densities, the grey-level function on a scintigram is heavily polluted by a 
quantum-type noise, where the size of the error in a particular point can be relatively 
large. Eliminating noise represents a major problem in the analysis of this type of 
data sets, and a standard approach in image analysis is to use a smoothing convolution 
filter before processing the images. Persistence cancelling provides a different way of 
avoiding this problem. Figure 18.41 shows, how the algorithm traces the maxima of the 
grey-level function on unfiltered images. The persistence level has been set to 50 (on 
a scale from to 255). 

Our experiment shows that running the algorithm on images smoothed by a Gauss- 
ian filter and setting the persistence to a lower value of p = 10 (an even lower value, 
for example p — 0, gives too many local maxima which correspond to outliers despite 
the smoothing, so it is not good for tracing the locations) produces a similar pattern 
as the one in Figure 18.4} but typically with only one peak within one kidney. Multi- 
ple peaks like those appearing at times t — 2, t — 8 and t = 9, are also of interest, 
though, since they give insight into the activity within the kidney. So prefiltering 
images seems to cause loss of useful information and using unpreprocessed images 
and higher persistence is a better choice from this point of view. 
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Figure 8.2. Tracing cavities through head CT scans 



Persistence also enables analysis of images at different scales. For example, once 
the kidneys are located in the above experiment, setting persistence to a lower level 
gives information on the distribution of the agent within the kidneys, and therefore 
information on the functioning of parts of a kidney. 

Finally, we note that the algorithm has also been used in a machine learning algo- 
rithm from a robotic domain where a robot observing a red and a blue ball with an 
on-board camera is taught the concept of occlusion from examples of image sequences 
jZBJPSj . 
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Figure 8.3. Simulated cavities in the head 




Figure 8.4. Tracing the maxima of radioactivity in scintigrams 
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